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Abstract 

We provide an algorithm to compute the 2-norm maximum of a multilinear map over a product of 
spheres. As a corollary we give a method to compute the first singular value of a linear map and 
an application to the theory of entangled states in quantum physics. Also, we give an application 
to find the closest rank-one tensor of a given one. 
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Introduction. 

A lot of problems in mathematics need to maximize a bilinear form over a product of spheres, 
for example the 2-norm of a matrix is given by the maximum of the bilinear form (x, y) — > x'Ay, 
where ||x|| = \\y\\ = 1. Another interesting problem is to find the closest rank-one tensor of a 
given tensor £ aijkXi®yj®Zk- To answer this problem one has to find the maximum of a trilinear 
form over a product of three spheres (see the examples). 

This article provides an algorithm to find the maximum of a multilinear map over a product 
of spheres, 

f:i" l+l x...xi"' + Ur tl+l , max \\€(x u . . . , x r )\\. 

I|a-iII=...=I|.v,.||=i 

We have reduced the problem of finding the maximum of I to a problem of finding fixed points of 
a map W : P" 1 x . . . x P" ,+1 — > P" 1 x . . . x P" r+1 . The advantage of this reduction is the possibility 
to count the number of extreme points of (, and also, to find the fixed points of W solving a 
system of polynomial equations. There are standard algebro-geometric tools to solve systems of 
polynomial equations. 

In lSection ll we review some concepts and definitions in algebraic geometry, such as, projec- 
tive space, maps, products of projective spaces and maps between them. We use these definitions 
in the article. 
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In lSection 21 using Lagrange's method of multipliers, see [Q1 §13.7], we reduce the problem 
of finding the maximum of a multilinear map I, to the problem of finding fixed points of a map 
W. We compare our approach with the ones in the literature. 

In lSection 3 1 we make a digression to discuss the number of extreme points of a multilinear 
map over a product of spheres. We use intersection theory to count the number of fixed points 
of the map W : P" 1 x . . . x P"' +1 P" 1 x . . . x P"' +I . Recall that the number of fixed points of 
a generic map F : P N — > V N of degree d is 1 + d + . . . + d N . In this section we give a formula 
to compute the number of extreme points of a multilinear map over a product of spheres. If the 
map is generic, this number is achieved over C, and if it is not generic, this number is a bound 
when the extreme points are finite. In the literature, the extreme points of £ are called singular 
vectors (see [ 16]) and in this section we count them. 

In lSection 41 we use our approach to find the maximum of a bilinear form over a product of 
spheres. In the bilinear case, the map W, induces a linear map L : P N — > P", where Nisa natural 
number, and we prove that for a generic q e V N , the sequence \q, L(q), L 2 (q), . . .} converges to 
the absolute maximum. In other words, the absolute maximum is an attractive fixed point of L. 
Also, with the same tools, we give an algorithm to find the spectral radius of a square matrix. 

In ISection 51 we use the theory developed to present the algorithm. We take advantage of 
a result in ISection 2l the classes of extreme points of the multilinear form I are in bijection 
with the fixed points of the map W. We reduce the problem of finding fixed points of W to 
solve a system of polynomial equations with finitely many solutions. In the literature about 
computational aspects of algebraic geometry, there exists a lot of algorithms to solve a system 
of polynomial equations with finitely many solutions, see [9]. This gives us the ability to find 
the absolute maximum of t. It is important to mention that the system of polynomial equations 
obtained with our approach is slightly different from the system of polynomial equations obtained 
naively from the method of Lagrange's multipliers. Our approach in projective geometry, allows 
us to find the correct solution removing some constrains. In the first part of the section, we 
present a direct method to find the maximum value of a generic multilinear form over a product 
of spheres. Basically, it reduces to find the spectral radius of a matrix. In the second part of the 
section, we give an algorithm to find the point (jci, . . . ,x r ) e R" 1+1 x . . . x R"' , where ||x,j| = 1, 
1 < i < r, such that \6(x\, . . . , x r )\ is maximum. This last algorithm, requires the extra hypothesis, 
2«i , . . . , 2n r <«! + ... + n r . 

In ISection 6l we use the theory developed to compute a lot of examples and applications. 
One of them is the ability to find the closest rank-one tensor of a given tensor. We prove that this 
problem is well posed and we apply our algorithm to solve it. Another application that we will 
give is related to quantum physics. It is a criterion of separability, given a quantum state, we can 
say if it is separable (see Remark 22l for definitions and related concepts). 



1. Review on Projective Geometry. 

In this section we give basics definitions that we are going to use, such as, projective space, 
maps, projective tangent space, product of projective spaces and maps between them. In this 
section we are assuming that the base field is R, but all the definitions are true in the complex 



case. All the notions in this section may be found in 1113 1 



Definition 1. Let n be a natural number and let R" +1 be a real vector space of dimension n + 1. 
The projective space, P", is the space of lines passing throw the origin in R" +1 . We say that the 
dimension of P" is n. Every nonzero vector v in R" +1 determines the line [v] that joins v with the 
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origin e R" +l . The vector v and Av, for A e R, A £ 0, determines the same point [v] e P". 

Let's fix a basis {vo, . . . , v„] of If the coordinates, in this basis, of v are («o, . . . , a„), 

then the coordinates of the point [v] are 

[v] = (ao : . . . : a„) = (i«o : ■ • ■ : ^««), A e R, 1 ^ 0. 

In general we denote [v] e P" to remark that the point [v] is represented by the vector v e K" +1 . 
Also, we denote an arbitrary point in projective space, p e P". The projective space P" is a 
compact space. 

Let n and m be two natural numbers. We say that a polynomial P in n + 1 variables is 
homogeneous of degree d, where <i is a natural number, if 

P(Ax , . . . , Ax n ) = A d P(x Q , . . . , x n ), A e R, /I # 0. 

For example, a linear form is homogenous of degree 1 . 

A map F from P" to P m , denoted F : P" — > P" 1 , is given by m + 1 homogeneous polynomials, 
F , . . . , F m , of degree <i 

F = (Fo : . . . : F m ) : P" — > P m , F(x) = (F (x) : . . . : F m (x)), x £ P". 

The homogeneity of the polynomials Fo, . . . , F m , implies that the value of F at [v] and at [Av] is 
the same in P m . We say that F has degree equal to c/. When d = 1 we say that F is linear. 

Let m, . . . ,n r be a list of natural numbers. A multihomogeneous polynomial is a polynomial 
F in variables x' Q , . . . , xj,,, for 1 < i < r, such that 

FOli* 1 , . . . , /U r ) = ^ 1 . . . AfP(x l x r ), x l = (x , . . . , 

The vector (Ji , . . . , c/ r ) is called the multidegree of F. For example, a multilinear form is a mul- 
tihomogeneous polynomial of multidegree (1, . . . , 1). 

A map F : P" 1 x . . . x P"' — > P m , where m e N, is given by m + 1 multihomogeneous 
polynomials, Fo, . . . , F m , of multidegree (di,..., d r ), 

F(x 1 ,...,x r ) = (F (x 1 ,...,x r ) : ... : F m (x\ . . . , x r )), x'' e P"', 1 < i < r. 

The multi-homogeneity of the polynomials Fo, . . . , F m , implies that the value of F at ([v 1 ], . . . , [v r ]) 
and at ([Aiv 1 ], . . . , [A r v r ]) is the same in P m . We say that F has multidegree (d\,..., d r ). 

Finally, a map F : P" 1 x . . . x P"'" — > P mi x . . . x P m - is given by s maps F = (F u . . . , F. s ), 

Fj : P" 1 x ... x P"'" — > P'"', 1<;<5. 

Note that the multidegree of F, may differs from the multidegree of Fj, i + j. When all the forms 
{Fi , . . . , F s ] are multilinear, we say that F is a multilinear map. 
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Definition 2. Let n and m be two natural numbers and fix a basis for M" +1 and for R"' +1 . Every 
vector v e has associated a vector space of dimension n + 1; the tangent space, denoted 

r v M" +1 . 

A polynomial map F = (Fo, . . . , F m ) : — > R'" +1 such that F(v) = w determines a linear 
map, dF v , called the differential of F at v, 



dF v : T v R n+1 r u ,R m+1 , rfF v (a , . • • , a n ) 



u=o 0X1 i=o 0X1 



In projective space the situation is similar, lfl3i p. 181]. Every point x eP" has associated an 
«-dimensional projective space; the projective tangent space, denoted T X P". A map F = (Fq : 
. . . : F,„) : P" — > P m of degree d such that F(x) = y induces a linear map between projective 
tangent spaces, 



dF x : T X V -> T,P m , dF x (a :...:«„) 



ZdF •A dF m 



V !=0 



t=0 



Given that the partial derivative of a homogeneous polynomial is also homogeneous, the map 
dF x is well defined. 

Remark 3. Recall the Euler relation for a homogeneous polynomial P of degree d, lfl3l p. 182], 



ZdP 
— (v)v; = d-F(v), v = Oo, . ..,vjv) e. 

*=o dxi 



The relation follows at once by differentiating both sides of the equation P(Av) = A d P(v). 

If F — (Fq : . . . : F#) : P w — > P w is a map of degree d and jc € P N is a point such that 
F(x) = x, then, using the Euler relation, we get dF x (x) = x, 



1 n 



dF x : T X P N -> T X P N , dF x (a : . . . : a N ) = 



\i=0 

In particular, if the vector v e M. N+1 represents x e P", x — [v], and the matrix dF x represent the 
linear map dF x , 

3 -p 

(dF x ) =^r 1 (v), 0<i,j<N, 

then, v is an eigenvector of dF x . Let's compute the eigenvalue of the eigenvector v. Given that 
F(x) = x there exists a nonzero real number A such that (Fo(v), . . . , Fjv(v)) = Av. Then 

1 N dF 

^ = F j{v) = -Y J ^liv)v i , 0<j<N. 

/=() 

Then, the eigenvalue of v is d ■ A, where d is the degree of the map F. 



dF 0, \ V dF N 

-i ox-i £a dxi 



i=0 
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2. Theory for a multilinear map. 



In this section we translate the problem of finding a maximum of a multilinear map to a 
problem of finding fixed points. Let's present the notation and some basics preliminaries. 
Let S" be the sphere in R" +I , 



S" = {u e R" +1 : \\u\\ = Vl«ol 2 + 777 + \u n \ 2 = 1), 

and let (— , -) : R" +1 x R" +1 — > K be the inner product, (x,y) = x^yo + . . . + x„y„. The norm 
associated to this inner product is the usual 2-norm, (u, u) = ||m|| 2 . 
When the codomain of a map is K, we say that the map is a form. 

Lemma 4. Given a multilinear map I : R" 1+1 x . . . x R" r+1 — > R J+1 there exists a multilinear 
form I, 

t : R" 1+1 x ...xR"' +1 xR' s+1 — > R, 2(x u ■ ■ ■ , x r ,y) = {t(x u . . . , x r ),y), 



such that 



max \\({xu . . . ,x r )\\ = max \t(xi, . . . ,x r ,y)\. 
||jtill=...=IWI=i ||.t,||=...=|U,-ll=llyll=i 



Proof. The proof is bases on the compactness of the sphere. Let (jti, . . . , x r ) e S" 1 x . . . x 8"' be 
a point such that z = €(x\, . . . , x r ) has the maximum norm and let y = z/\\z\\- Then 

\£(x u . . .,x r ,y)\ = |<z,y>| = = \\z\\ = U(x u ■ ■ .,x r )\\ 

\\z\\ 

max \i(xi, . . . ,x r , y)\ > max ||£(jci,...,jc r )||. 

l|A',ll=...=lk,.||=||v||=l Hjc,||=...=||i,||=1 

Analogously, let (x\,.. .,x r ,y) e x . . . x x § s be a point such that \({x\, . . .,x r ,y)\ is 
maximum. Let z = t(x\, . . . ,x r ). Then, 

\i(x u . . .,x r ,y)\ = |< z ,y>| < IklUMI = \\£(x u . ■ .,x r )\\ 



max \{(xu ■ ■ ■ ,x r ,y)\ < max \\£(xi, 

II-V1||=... = ||.V,.||=||V|| = 1 ||. tl ||=... = ||A-,.||=l 



□ 



As a corollary of the previous lemma, we will work with multilinear forms. Specifically, to 
make the notation easiest, we will work with Z : R" +1 x R m+1 x R I+1 — > R a trilinear form. Our 
goal is to find the maximum of t over a product of three spheres. 

Using Lagrange's method of multipliers, ([Ql §13.7]), we know that the extreme points of t, 
over §" x x S s , satisfy 

( d£/dxi(x ,...,x n ,y Q ,...,y m ,ZQ,...,Zs) = 2ax ( , < i < «, 
| de/dy j(x Q , ...,x n ,y ,...,y m ,Zo,---,Zs) = 2/3yj, < j < m, 
[ d€/dzk(xo,...,x„,y ,...,y m ,zo,---,Zs) = 2Az k , < k < s, 

a,{3,AeR, \\x\\ = \\y\\ = ||z|| = 1. 
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Let's use a better notation, 



x = (x ,...,x n ), y = (y ,...,y m ), z = (zo,...,z s ), 

— (x,y,z) = — (x,y,z),..., — (x,y,z) ], ^-(x,y,z) = —(x,y,z),...,- — (x,y,z) , 
ox \dxQ dx„ j dy \dy dy m j 

dt, , [dt , dt , \ „„, [dt , dt, dt, \ 

— (x,y,z) = —{x,y,z),...,—(x,y,z)\, \t(x,y,z) = —(x,y,z),—(x,y,z), — (x,y,z) ]■ 
oz \ozo ozs I \ox oy oz I 

Definition 5. A point (x, y, z) e S" X § m x § ! is called an extreme point of I if it satisfies the 
system of equations 

Vt{x,y,z) = (2ax, 2/3y, 2Az), 

for some a,/3,Ae M. Note that if (x, y, z) is an extreme point, then (+x, ±y, +z) is also an extreme 
point. We say that they belong to the same class. 

Proposition 6. There is a bijection between classes of extreme points oft and fixed points of the 
map 

W : F" x P'" x P' 5 -> P" x P m x P' ! , 
(W, M, M) ^ f?(W, M, [z]). ^(W, [y], [z]), ^(W, [y], [z]) 

\ ox oy oz 

Proof. Given an extreme point (x,y,z), consider ([jt], [y], [z])- This assignment is independent 
of the class of (x, y, z). By definition, it gives a fixed point of W. 

Given a fixed point ([x], [y], [z]) e P" x P m x P 1 of Vt, consider representatives x,y,z such 
that ||x|| = \\y\\ = \\z\\ = 1. Then (x,y,z) is an extreme point of t. □ 

Remark 7. The map V/" : P" x P m x P J -> P" x P'" x P' ? from |Proposition 6| is not defined over 
the closed subset 

dt dt dt 1 

(x,y,z) I — (x,y,z) = or — (x,y,z) = or — (jc,y,z) = \ c P" x P" x P s , 
ox oy oz J 

but this set is empty if and only if the hype rdete rminant of € is zero. The hyperdeterminant is a 
polynomial in the coefficient of t, for the definition and some properties see II12L §14]. 
By a result in HI §14, 1.3], if 

2n, 2m, 2s < n + m + s 

then a generic choice of t will make W defined everywhere. 

In the article fltSll . there is a definition of singular values and singular vectors for a multilinear 
form t. For example, for a trilinear form t, the author defined the singular vectors of t as the 
solutions of the system Vt(x,y,z) = (2ax,2/3y,2Az). It is the same as our definition of extreme 
points. It is of interest to know the number of singular values/vectors of t, and in lSectiorTJl we 
count them. In the same article, the author proved that the first singular value is the maximum of 
t over a product of spheres. Also, under the hypothesis 2n, 2m, 2s < n + m + s, he proved that the 
hyperdeterminant of t is zero if and only if is a singular value of t. The hyperdeterminant is a 
polynomial in the coefficients of t, so, in the generic case, the number is not a singular value 
oft. 

There exists another article to mention, 17[]. In it, the authors proposed a multidimensional sin- 
gular value decomposition, but it does not preserve the properties that we need, for example, that 
the first singular value of t corresponds to the maximum of t over a product of spheres. 
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3. Number of extreme points of a multilinear form. 

In this section we use Intersection Theory ( fill 8.4]) to count the number of fixed points of 



a generic map P" x P m x P s —> P" x P m x P s over C. Recall from Proposition 6 that there is a 
bijection between fixed points of 

W : P" x P'" x P' s — > P" x P" x P s 

and classes of extreme points of the trilinear form £ over S" x S m x S ! . The number of fixed 
points of W gives a bound to the number of classes of extreme points of £ that contains a point 
with maximum value. It is known that if F : P N — > ¥ N is a generic map of degree d, then F has 
1 + d 2 + . . . + d N fixed points, ifToL 1.3]. Here we generalize this result to a generic map between 
products of projective spaces. 

Before we continue with this section, let's make a survey of some related concepts that are in 
the literature. 

In |0], 1 21], l20ll . 01, H5I1 and Jst] there is a notion of eigenvectors and eigenvalues associ- 



ated to a multilinear form i. There are a lot of applications and in Q5[| , the authors counted the 
number of eigenvalues of I as the number of roots of a characteristic polynomial associated to 
I. The idea is to look at I : C" x . . . x C" C as a polynomial map P : C" (C") v = C", 
P(x) = l(x, . . . , x, — ) and then, an eigenvector of I is a vector x e C" such that P(x) = Ax. If ( is 
m-multilinear, P has degree m - 1 and as a map P"~ 1 — > P"~ 1 it has (m - 1 )"~ 1 + (m - 1 )"~ 2 + . . . + 1 
fixed points, i.e. eigenvectors of P. They arrived at this number using toric varieties and Newton 
poly topes. 

In il 7.1.4] and lfl9[ 3.1] there is a theory of multihomogeneous Bezout number, or m- 
Bezout. The m-Bezout gives an upper bound on the cardinality of the intersection of multiho- 
mogeneous polynomials in P" 1 x . . . x P"' . Given that we are counting the fixed points of a map 
F : P" 1 x . . . x P" 1 — > P" 1 x ... x P" A , in order to apply this formula, we need to realize the fixed 
points of F as an intersection in some product of projective spaces. Concretely, the intersection 
of the graph of F and the diagonal. Let's make an explicit example. Assume for simplicity, that 
F is linear, F : P" — > P", we will see that the m-Bezout formula gives a very bad bound. Recall 
that the number of fixed points in this situation is the number of eigenvectors, that is, n + 1. Let's 
apply the formula to the equations of the graph T = {(x, F (x))} and the diagonal A = {(x, x)}. The 
points in the intersection satisfy the following equations, 

((x Q x n \ (y : . . . : y n )) e P" x P", yi Fj(x) = yjFfc), x iyj = Xj y u < i, j < n. 

Note that the equations correspond to the fact that the following matrices have rank one, 

yo ■■■ y n \ 1*0 ■•■ x„ 

F (x) ... F„(x)J' [yo ... y„ 
By abuse of notation, we denote the equations, 

(x,y)€P"xP", y = F(x), x = y. 

Given that the equations have bidegree (1,1), the m-Bezout number is the coefficient of a^a^ 1 
in the polynomial (a\ + a2) 2n+2 . It is the binomial ( 2 ^' 4 + 1 2 ) +n + \. 



Bernstein proved in [4] that the number of solutions of a sparse system equals the mixed 
volume of the corresponding Newton polytopes. A sparse system is a collection of Laurent poly- 
nomials, 

fi= Ci .vi.-,v„*? ■ • ■ x n > i=h---,n 

(Vi,...,V„)£^i 

where Mi are fixed finite subsets of Z". Its convex hull Q t = conv(J?l,-) c R" is called the Newton 
poly tope of fj. Consider the function 

R(A U . ..,A n ) := vol(Aifii + . . . + A n Q n ), A, ; > 0, i = 1, . . .,n. 

where vol is the usual Euclidean volume in R" and Q + Q' denotes the Minkowski sum of poly- 
topes. It is a fact that R is a homogeneous polynomial and the coefficient of the monomial 
/li... A„ is called the mixed volume of Q\, . . ., Q„. The mixed volume (i.e the number of solu- 
tions of a sparse system) is a very difficult number to compute by hand. In some situations, this 
is possible and in the general case, there are a lot of algorithms to compute it. In our situation, 
we are working with a multihomogeneous polynomial system, and using Bernstein's theorem, in 



the paper [18], the author gives a recursive formula to compute this number. In fact, it is proved 
that, under some hypothesis, if the system is over R and the functions are generic, then all the 
solutions are reals. Here, we present a different and more direct method using intersection theory. 



In the article 11161 3] there is a definition of generalized singular values for a generic multilin- 
ear form. In this section we count them. 



Let's make an introduction to Intersection Theory. The germ of intersection theory is the 
Fundamental Theorem of Algebra. It implies that given a generic homogeneous polynomial in 
two variables F of degree d, the set of zeroes \x e P 1 , F(x) = 0} has d points. Generalizing this 
idea, Bezout's theorem, says that given two generic homogeneous polynomials in three variables 
of degree d and e, the set of zeroes \x e P 2 , F\{x) - F2(x) = 0} consists of de points. In P r 
the situation is similar, if F\ , . . . , F r are generic homogeneous polynomials of degree d\,...,d r 
respectively, the intersection has d\d% . . . d r points. 

To formalize these ideas, let's introduce the Chow ring of P r , ifTll proof of Prop. 8.4] 

A(P r ) = Z[cr]/(c/ +I ). 

Every variety XcF has a class, [X] e A(P r ). The intersection of two generic varieties X C\Y 
corresponds to the product of the classes [X].[Y] = [X n y]. Two different varieties may corre- 
spond to the same class, for example, every hypersurface of degree d corresponds to the same 
class, da, where a is the class of a hyperplane. For example, a r corresponds to the intersection 
of r generic hyperplanes, i.e. a point. The product 

(d\a).(d2a) (d r a) — d\ . . .d r a r 

corresponds to the intersection of r generic hypersurface of degree d\,...,d r respectively. We 
get d\ . . .d r points in the intersection as mentioned. The class of a variety of codimension c is a 
homogeneous polynomial of degree c in Z[a]/{a r+l ). 

The Chow ring is very useful to solve problems in enumerative geometry. For example, to 
count the number of fixed points of a generic map P' — > P'' the procedure is the following. Let 
A(P r x P') be the Chow ring of P r x P r , it is A(P r x P') = Z[a, a]/(a r+ \ a r+1 ), H2 Ex. 8.4.2]. 
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Let [A] e A(P r X P r ) be the class of the diagonal, A = {(*, x)}, and let [r] € A(P r x P r ) be the 
class of the graph of F, T = {(x, F(x))}. Given that 

dim A + dimT = dim(P r x P r ), 



the product [A].[T] is a multiple of the class of a point, da r a r , 111 ll §8.3]. The coefficient d is the 
number of fixed points of F. 

The Chow ring of a product of projective spaces, Ull Ex. 8.3.7], is 



A(P ni x...xP') = A(P'")< 



■A(F n *) = Z[a l ,...,a k ]/(a> 



Note that in A(P"' x . . . x P"*) there is only one class of a point, a" 1 . . . a" k k , so there is a well 
defined map called degree. The degree of a class is the coefficient of a" 1 . . . a" k k . It may be zero 
(or negative). 

The last thing to mention is that every map F : P" 1 x . . . x P" A -» P"" x . . . x P"" induces a 
morphism of rings, 111 ll Prop. 8.3 (a)], 



F* : AQT 1 x . 



') — > A(P"' x ... x P" 1 ), F*([X]) = [F~ 1 (X)]. 



For a more extensive treatment of intersection theory, see [ 14, §A], [11]. 



Let's use the previous introduction. First, we will compute the number of fixed points of 
a generic map F : P r — > P r of degree d. Then, we will adapt the proof to a generic map 
F : P"xP'"xP s — > P"xP m xP 5 formed by multihomogeneous polynomials of some multidegree. 

Proposition 8. The number of fixed points of a generic map F : P r — > P r of degree d is 

l+d+... + d r . 

Proof. The following proof is standard in intersection theory. The fixed points of a map F : 
pr pr ma y k e computed in A(P' x P' ) as the degree of the product of the class of the graph of 
F, [T], and the class of the diagonal, [A]. First, let's find out the class of the diagonal, 

[A] € A'(P r x P r ) = Z[a, a]/(a r+ \ a r+l ). 

Being of codimension r, the class is a homogeneous polynomial of degree r, 

[A] = t {) a'~ + t\aa r ~ l + . . . + t r -\a r ~ x a + t r a r , f ; e Z 

Here, a represents a class of a hyperplane in F r and a' represents the intersection of i of these 
generic hyperplanes, in other words, a' is a generic space P r ~' inside P r . Same for a and aK 
Viewed in P r x F', a' is the class of U x V, a' = [U x P r ], and a* is the class of W x V, 
a> = [F r x V], where dim U = r — i and dim V = r — j. The class da' , represents a product of 
general linear spaces U x V c P r x P r , where dim U = r — i and dim V — r - j. 

The class of the diagonal is determined by the coefficients fo, . . . , t r . Note that f, = [A].a r ~'a l . 
Then, we need to count the number of points in (U x V) n A, 

( U x V) n a s u n v = {p} 
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t = . . . = t r = 1 => [A] = V a i a r '. 

!=0 

Now, let's compute the class of the graph of a map, F = {(x, F (x))} cf'x P r , 

[H e A r (P r x P r ) = Z[a, a]/(a r+1 , </ +1 ), 

it is also a homogeneous polynomial of degree r, 

[r] = Toa r + T\aa r ~ l + . . . + T r -\a r ~ x a + T r a r , t ; e Z. 

Again, we have t, = [T].a r ~'a' so we need to count the points in T n (U x V), where dim t/ = i 
and dim V = r-i, 

rn(UxV) = {xe U\F(x) e V) = U n F _1 (V) c P r . 

If F is formed by homogeneous polynomials of degree c/, the pull-back of a hyperplane is a 
hypersurface of degree d, then 

[U n F _1 (V)] = o i .F*(tO = a , '.F*(Qr) r - i = aXda)^ = *TV. 

Then, t/ n F~\V) has d'"' points, i.e. t,- = cT ! ', 

[r] = + d^W" 1 + . . . + da^a + a r => 

r r r r 

[A].[r] = fl'o^'Xj] d r - j a i a r - j ) = ^ rf r -V' + V r - (!+ ;' ) = J] ^ = 1 + d + . . . + d r . 

Given that a constant map has one fixed point, we use the convention d () = 1 for d = 0. □ 

Let's adapt the previous calculation to P" xP m xP*. We will compute the class of the diagonal 
and the class of a graph, and then we will multiply them to obtain the number of fixed points. 

Theorem 9. The number of fixed points of a map F = (F U F 2 ,F 3 ) : P" x P m x P s -» P" x P m x P' ? 
is the coefficient ofa"j3 m y s in the following polynomial in Z[a,f), y], 

n m s 

X Z Yj {dia + W + ^rrX^ia + e 2j 6 + e 3 y) m -j(fa + fifi + fiyY^aW, 

i=0 7=0 k=() 

where {d\,d2,d 3 ), {e\, ei,e{) and (fi,fi,f 3 ) is the multidegree ofF\, F2 and F 3 respectively. 
For a generic map F : P" 1 x...xF" — > P" 1 x...xP" the result is similar. 

Proof. The class of the diagonal A = {(x,y,z,x,y,z)} e P" x P m x P' ? xPx P ffl x P s , is a 
homogenous polynomial of degree n + m + s, 

[A] e A(P" x P m x P s x P" x P m x P' ? ) = Z[a,/3,y,a,b,c]/(a" + \/3 m+1 ,y s+ \a n+1 ,b m+ \c s+1 ), 

Instead of doing the same computation as before, let 

tt ia : P" x P m x P s x P" x P m x P s -> P" x P", 
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be the projection in the first and fourth factor (same for 7:2,5 and ^3^) and let A„ cP'x P" be the 
diagonal of P" (same for A,„ and A s ). Then we have 



n m s 



[A] = 7r* 3 ([A„]).7r* 5 ([A m ]).7r 3 * 6 ([A s ]) = £ J] £ ct^VfT^c^. 

i=0 j=0 k=0 

The class of T = {(x,y,z,F(x,y,z))} Q P" x F" x P s x P" x P'" x P'\ is a homogeneous 
polynomial of degree n + m + s, 

[T]= J] TifttfHaWVprty, T fJW/t ,eZ=> 

i+ j+k+i'+ j'+k'-n+m+s 

n m s 

deg([A].|T]) = £££Ty^. 

i=0 7=0 it=0 

Where t/eg is the coefficient of a"a"b"'f3 m c s y s . Note that the integer t,^,-^ is the degree of 
[r].a n ~'a l W'-ipj c s ~ k y k ; the number of points in 

r n (t/i x t/ 2 x t/ 3 x Vi x v 2 x y 3 ), 
c/i, Vi c p", u 2 , v 2 c p m , c/ 3 , y 3 c p s , 

dim {/ 1 + dim V*i = n, dim f/2 + dim V2 = m, dim {73 + dim V3 = s ==> 
r n (Ui x C7 2 x C7 3 x V\ x V 2 x V 3 ) =* (f/i x C7 2 x C7 3 ) n F _1 (Vi x V 2 x V 3 ) c P" x P m x P s . 
Let's use the fact that F is equal to (Fi,F 2 , F3), 

F] : P" x P'" x P' ? -> P", F 2 : P" x P" ! x P' ? -> ¥"', F 3 : V x ¥"' x P s -> P" 

where (d\,d 2 ,di), (e\, ei,ei) and if\,f 2 ,fo) is the multidegree of F\,F 2 and F3 respectively. 
Then 

F~ l {V, x V 2 x V 3 ) = F^CVi) n F- 2 \V 2 ) n F 3 - 1 (V 3 ). 
Thus, the class of the intersection that defines t,^,^ in the Chow ring A(P" x P" 1 x P s ), is 

atp-fidia + ^ + *r)"" ! '(eia + e 2 /3 + e 3 y)" W (/i a + / 2 j6 + /sT)*"* 

□ 

Example 10. Let's apply the previous formula to W where I : § 2 x § 2 x § 2 — > R is a generic 
trilinear form. The multidegree of dl/dx, di/dy and cV/<9z is (0, 1, 1), (1,0, 1) and (1, 1,0) re- 
spectively. Then the number of fixed points of this map (over C) is equal to 37. The number 37, 
according to Ilia, 3], is the number of generalized singular values of I. 

Example 11. Let's apply the formula to count the number of eigenvectors of a generic linear 
map L : R" +I -> R" +1 . The map L induces a map P" -> P" of degree 1, then 



Y = V or" = (n + l)a". 



The map P" — > P" has n + 1 fixed points over C, that is, L has « + 1 eigenvectors over C. 
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Example 12. Finally, let's apply the formula to find the number of singular values of a generic 
linear map L : R" +1 -> M" ,+1 . The map L induces a bilinear form t : M" +1 x R'" +1 -> K and the 
bidegree of d£/dx and dt/dy is (0, 1) and (1,0) respectively (assume n > m). 

n m n m m 

YjYu^^^F = Yj Yj 13 ' 1 ^"" 1 ^ = Yj/3"- ( "- m+i)+i a m - jHn - m+j) = (m + l)a"j8 m . 

i=0 j=0 i=0 j=0 7=0 

Then, W : P" x P" -> P" x P m has m + 1 fixed points, that is, L has m + 1 singular values over 
C. We used the variational definition of singular values, see flol . In case n < m we can use the 
fact that the number of non-zero singular values of L and V : — > PJ ,+1 are the same. 



4. Theory for a bilinear form. 

In this section we present a method to find the maximum of a bilinear form, t , over a product 
of spheres, S" x § m . This case is very special and the method presented here does not work for a 
general multilinear form. 

The key point of this method is the fact that the partial derivatives of i — 2 a^x/yj are linear, 

— (x,y) = £(ei,y), — (x,y) = e(x,ef), (x,y) e R" +1 x R m+1 , < / < n, < j < m, 
ox-i dyj 

where e,- and ej are canonical basis vectors of K" +I and R m+1 respectively. The map Vt induces 
a linear map L : P"+'» +1 pB+m+i_ Let (x : . . . : ;t;„ : y : . . . : y m ) be a point in F n+m+1 . Then 

L(x : . . . : x n : y : . . . : y m ) = (t(e ,y) : t(e n ,y) : £{x, e ) : £{x, e m )) . 

This map is well-defined. Let A e M, A # 0, 

L(Axo : ... : Ax n : Ay : ... : Ay m ) = (£(e Q , Ay): t{e„, Ay) : {{Ax, e ) : . . . : ({Ax, e m )) = 

{A{{e ,y) : At{e n ,y) : Ai{x, e ): A({x, e m )) = L(x : . . . : x„ : y : . . . : y m ). 

Theorem 13. Let p — {x,y) eS"x S m be an absolute maximum of €. Then 

lim U{q) = [p] 



for a generic q e 



Proof. Let A e R«+ m +2xn+m+2 ^ e a jna^-jx representing the linear map L. Given that L is linear, 
the differential of L at any point, q, is equal to L, 

dL q = L, V ? eP" +m+1 . 

In particular, the matrix A, also represents the differential of L at p, 

A = dL p . 

Let {vo, . . . , v„ +m+ i) be a basis of H> n+m+2 formed by eigenvectors of A. Let Aj be the eigenvalue 
of Vj, < i < n + m + 1 . By IRemark 31 we know that p is an eigenvector of A with eigenvalue 
t{p). In particular, if the magnitude of Aq is maximum, then \ Aq\ — \({p)\ and [p] = [vq]. 
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Let z G R' !+m + 2 be a vector representing q such that z = flo^o + • • • + tfn+m+iV n+m+ i, ao + 0. 
Then 

n+m+l n+m+\ n+m+l 

L r (q) = [A r . ^ djVi] = [a A r Q v + ^ fl,/l[vj = [a vo + 2j fl <'7" v <] ~* [ a « v o] = P- 
i=o i=i i=i 

□ 



In the proof of the previous theorem, we saw that the iterations of a linear map in projective 
space converges to an eigenvector of maximum eigenvalue. In particular, given a square matrix 
A e p"+ lx " +1 and a generic vector w e R" +1 , the sequence {[w], [Aw], [A 2 w] . . .} c P", converges 
to a point [v]. The vector v satisfies Av = Av, where \A\ is spectral radius of A. 

The rate of convergence of this method is linear. 

Remark 14. Let's give an algorithm to find the absolute maximum of a generic bilinear form, 

t : R" +1 x R" ,+1 — > R, 

based on ITheorem 131 Let W = (d£/dx,d£/dy) be the gradient of £ and let q be the initial 
condition, where q = (x,y), x e R" +1 , x + and y e R" !+1 , y + 0. 

Input: A bilinear form I : R" +I X M m+I R. 

Output: The absolute maximum (a; )•) e S" X S m . 

1. Let q = q/\\q\\ and a!U = (1,0, . . . ,0). 

2. While a«A-)| is different from 1 do 

2.1 mix = q 

2.2 q = V£(q) 

2.3 9 = q/\\q\\ 

3. Let A : = (q , . . . ,q„), y = (q„ +i , . . . ,q„ +m+ i). 

4. Return {x/\\x\\, y/\\y\\). 

The iterations stops when the points in projective space are equal, in other words, when the 
cosine of the angle between q and aux is 1 or -1 (when they are aligned). Given that the absolute 
maximum is attractive (see lTheorem 131 ). the program ends. The maximum value is |^(;c,y)|. 

Remark 15. We may adapt this algorithm to a multilinear form, but in the multilinear case, in 
general, the absolute maximum is not an attractive fixed point. For example, the trilinear form 

£ : R 2 x R 2 x R 2 -> R, 

t(x,y,z) - 6x0^020 + IxxyoZo - faayiZo + 16x^^0 - \Ax y Q zi - \5x x y Q zi - llXoym + ^yizu 

induces a map P 5 — > P 5 of degree 2 without attractive fixed points. Even more, the 4-multilinear 
form f : R 2 x R 2 x R 2 x R 2 — > R, 

t(x,y,z, v) = Ax a y {) z v {) + 6x 1 >'oZ vo + x y,ZoVo - 6x 1 j z 1 v v - 5x y ZiV + 
Txiyizo - 5x yiZiV + 2x y z vi - 3xiy z Vi - 7x yiZ V| + 
9x l y [ zoVi - 9x y ziv l - 9xiy ZiVi - 6x y l ziv i + 8xi)>i£iVi, 

induces a map P 7 — > P 7 of degree 3 with two attractive fixed points. One is the absolute maxi- 
mum. 
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5. Presentation of the general algorithm. 



In this section we present an algorithm to find the maximum of a multilinear form over a 
product of spheres. First, we reduce the problem to a system of multilinear equations and then 
we resolve the system using algebro-geometric tools. In the first part of the section, we present 
an algorithm to find the absolute maximum of a multilinear form. In the second, we give an 
algorithm to find the point where the maximum occurs. This last algorithm requires some extra 
hypothesis. 

Proposition 16. Let t : R" +1 x R m+1 x R s+1 — > K be a generic trilinear form. There exists a 
bijection between classes of extreme points of I and solutions of the following system of trilinear 
equations in P" X F" X W, 

I£(xjei - x-,ej,y, z) = 0, < i < j < n, 
£(x, y jei - y,e j,z) = 0, <i < j < m, 
£(x, y, Zje-i - nej) = 0, < i < j < s, 

The vector e k satisfies (ek)i -Qifl + k and (ek)k = 1. 

In the multilinear case, we obtain a similar result; a system of multilinear equations. 



Proof. From Proposition 6 we know that every class of an extreme point of £ , is a fixed point of 
W : P" x F" x P s -> P" x F" x P s . If I is a generic trilinear form, we know that the number of 
fixed points is finite (see lSection 3i . 

A fixed point of W, ([jc], [y], [z]), satisfies 

d£/dxi(x,y,z) = 2axi, < i < n, 



{d£/dxj(x,y,z) = 2axi, < i < n, 
d£/d yj (x,y,z) = 2fyj, 0<j<m, 
d£/dz k (x,y,z) = 2A Zk , < k < s, 



where a, ft and A are three nonzero real numbers. In P" x P'" x P s , the equations are 

( xjd£/dxi(x,y,z) — Xid£/dxj(x,y,z), 0<i< j < n, 
< yjd£ldyi{x,y,z) = yid£/dyj(x,y,z), </<;'< m, 
{ Zjd(ldzi(x,y,z) = Zid£/dzj(x,y,z), < i < j < s, 

The result follows from the equalities, 

d£/dxi(x,y,z) = £(ei,y,z), d£/dyj(x,y,z) = £(x,ej,z), d£/dzk(x,y,z) = £(x,y,e k ). 



□ 



Let's present the algorithm to find the absolute maximum of a generic multilinear form. The 
resolution of the system is bases on Eigenvalue Theorem. Let's recall it. Consider a system of 
polynomial equations with finitely many solutions in C", 

fi(x\,...,x n ) = 

fm(xi, ■ ■ ■ , X„) = 
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where f\,..., fi„ are polynomials in C[x\ ,..,,x n ]. The quotient ring, 

M = C[x 1 ,...,x n ]/(f u ...,f m ), 

is a finite-dimensional vector space, Theorem 2.1.2]. The dimension of J{ is the number of 
solutions of the system. 

Every polynomial / e C[x\, . . . , x n ], determines a linear map M : — > M, 

M(g)=fg, g e C[x u ...,x n ], 

where g denotes the class of the polynomial g in the quotient ring M. The matrix of M is called 
the multiplication matrix associated to the polynomial /. 

Theorem (Eigenvalue Theorem). The eigenvalues of M are \f{p\), . . - ,f{p r )}, where [pi, . . -,pA 
are the solutions of the system of polynomial equations. See Theorem 2.1.4] for a proof. 



The algorithm in Appendix A first generates the following system of polynomial equations, 



xjd£/dxi(x,y,z) - Xid£/dxj(x,y,z), 0</< j < n, 

yjd{/dyi(x,y,z) = yid£/dyj(x,y,z), </<;'< m, 

Zjd£/dzi(x,y,z) = Zid£/dzj(x,y,z), < i < j < s, 

x 2 + . .. + x\ = 1, y\ + . .. + y 2 m = 1, zl + ... + z] = 1. 

Then, computes the real eigenvalues, {Aq, . . ., A r }, of the multiplication matrix associated to £. 
Finally, it returns A, such that > \Aj\ for all < j < r. This number, is the maximum of I over 
§" x x § s . 



For the algorithm and an implementation in Maple, see Appendix A 



Now, let's give an algorithm to find the point (x,y,z) 6 §" x §'" x § s such that \({x,y, z)\ is 
maximum. We need to use the following result (same notation as Eigenvalue Theorem), 

Theorem. Let x — A\X\ + . . . + A„x„ be a generic linear form and let M be its multiplication 
matrix. Assume that B — [I, Xu . . . , X„, . , .} is a finite basis of JL Then, the eigenvectors of M 
determine solutions of the system of polynomial equations. Specifically, ifv — (vq, . . . , v„, . . .) is 
an eigenvector of M such that Vo = 1, then (y\, . . . , v n ) is a solution of the system of polynomial 
equations. Even more, every solution is of this form. See 4£ §2.7.3/ for a proof. 

Note that the Theorem requires that the variables \x\, . . . , x„] are elements of the basis B. It 
could be the case that some variables are missing from B. For example, if X\, . . ., Xi e B, and 
jc<+i, . . . , x n i B, then, every missing variable, say Xj, is a linear combination of [x\, . . ., x,}, 

xj — aj\X\ + . . . + ajjXi, i + 1 < j < n 

If v = (1, vi, V2 . . .) is an eigenvector of M, the y'-coordinate of the solution corresponding to v, is 
Ojivi + . . . + ajjVj. See ||9l §2.1.3]. 

In order to apply the previous Theorem, we need to guarantee that the basis B contains all the 
variables. The affine system in Proposition 16 is, 



xjd£/dxi(x,y,z) = Xjd£/dxj(x,y,z), 0</< j < n, 
yjd£/dyi(x,y,z) = yid£/dyj(x,y,z), </<;'< m, 
Zjd£/dzi(x,y,z) = Zid£/dzj(x,y,z), < i < j < s, 
x = 1, yo = 1, Zo = 1. 
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The solutions of this system determine classes of extreme points of t . The genericity of i implies 
that all the extreme points of i, (x, y, z), satisfy xq + 0, yo + 0, zo + 0. Then, all the classes of 
extreme points appear as the solutions of the affine system in |Proposition 16| 

Theorem 17. Assume that I : x M m+1 x R s+1 — > E is a generic trilinear form and that 
2n,2m,2s < n + m + s, see \Remark 71 Then, the basis B of the affine system in Proposition 16 
contains all the variables. 

In the multilinear case, we obtain a similar result. 



Proof. Given that the equations in Proposition 16 are multilinear, the quotient ring, ft, is multi- 
graded. Let's denote ft{d 1 ,d i ,d i ) the multidegree part (d\,d2,d^), where d\,d2,di > 0. The hy- 
pothesis 2n, 2m, 2s < n + m + s, implies that the following set is empty, 

QP QP qp \ 

(x,y,z) e V" x W x P* | —(x,y,z) = or —(x,y,z) = or —(x,y,z) = (H = 0. 

ox oy oz J 

Then, the equations {d£/dxi}" =0 are linearly independent. Same for {di/dyj}'J =0 and {d€ I dzk\ s k _ - 
In the quotient ring, M, the partial derivatives, are proportional to the variables, thus, the variables 
are linearly independent too. For example, a basis for the multidegree part (0, 0, 0) is {1}, and a 
basis for the multidegree part (1,0, 0) is {xq, . . ., x n ). Even more, a basis for 

•^(1,0,0) © ^(0,1,0) © -^(0,0,1) 

is {x , . . . , x n ,y , . . .,y m , zo, . . . , z s }. 

Let's add the equations xq = yo = zo = 1 to the system of polynomial equations. The 
equations are not multilinear, so the corresponding quotient ring is not multi-graded, 

Jfl = J?l/<*o-l,yo-l,Zo-l>. 

Let's see that the variables {x\, . . . ,x„,y\, . . . ,y m ,Z\, . . . ,z s } are linearly independent in M. 
This implies that the basis B of 3\, formed by monomials, contains all the variables. 

n m s 

^ atXi + ^Bjyj + ^ A k z k = e M, a h Bj, A k 6 C => 
i=i j=i k=i 

n m s 

^ a iXi + J^Bjyj + ^ A kZk = P 6 ft, P e (x - l,y - l,Zo - 1), 

i=\ 7=1 k=l 

where P is a polynomial combination of xq - 1, yo - 1 and zo - 1. 

Denote P{d i ,d i ,ck) me multidegree part, {d\,d2, di), of P. Given that ft is multi-graded, we get 
the following equalities in ft, 

n m s 

P = P(l,0,0) + P(0,1,0) + ^(0,0,1)' 2u a ' X ' ~ °' 0) ' y^j^jyi = ^(0,1,0), / \hZk = P (0,0,1)' 

i=l 7=1 k=\ 

Using the fact that the variables \xo, . . . , x„) are linearly independent in ft, we obtain that xq is 
not a variable in P. Same foryo and zo- Given that P is a polynomial combination of xq - 1, yo - 1 
and zo _ 1, it must be 0. Then, a\ = . . . = a n = 0, B\ = ...= B m = and X\ = . . . = A s = 0. □ 
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Remark 18. The algorithm to find the point (x,y,z) e S" X x S s such that \i(x,y,z)\ is maxi- 
mum, is the following. 

The reader may adapt the algorithm to a multilinear form. 

Input: A generic trilinear form t : M." +i X K" ,+ I X M ,+1 -» M, 
where 2n + 2m + 2s < n + m + s. 
Output: The absolute maximum (x, y, z) € S" X S m X S s . 

1. Compute the system of trilinear equations of [Proposition 16] 

2. Add the equations xq = yo = Zo = 1. 

3. Compute a Gbbner basis for the resulting system, I. 

4. Find a basis B of <C[xi, . . . ,x„,yi, . . . ,y m ,zi, . . . ,Zj]/J. 

5. Compute the multiplication matrix of xi, 

6. Compute the eigenvectors of the multiplication matrix. 

7. For each eigenvector v do 

7.1 Normalize v such that v = (l,vi,...). 

7.2 Let x = (xi, . . . ,x n ) be such that .v, = Vm where 
cr, is the coordinate of .v, in B, 1 < i < n. 

7.3 Same for y and z- 

7.4 Normalize the points, x = Jt/||jt||, y = y/||y||, Z = z/l|z||. 

7.5 Evaluate fat (at, y, z) if the coordinates are real. 

7.6 Save the maximum. 

8. Return the maximum, (x, y, z)- 

In Step 5 of the algorithm we used the linear form X\ as a generic linear form. This fact is 
not restrictive. Given that the trilinear form is generic, we may suppose that the first coordinates 
of the classes of extreme points of I are all different. In other words, the eigenvalues of the 
multiplication matrix of X\ have multiplicity one. See §2.1.3]. 



6. Applications and examples. 



Let's start we some applications. First, we give applications of the iterative algorithm to find 
the maximum of a bilinear form. Then, we give applications of the general algorithm. 

Remark 19. Given a real matrix A, its first singular value (the 2-norm) is given by 

max x'Ay. 
IMI=lb>ll=i 

In other words, it is the maximum over a product of spheres of the bilinear form (x, y) — > x'Ay. 

An interesting fact of lTheorem 131 is that we can find the first singular vectors and the first 
singular value, |^(x,y)|, of I without using the spectral radius formula. Recall that the 2-norm of 
a matrix A is computed using the spectral radius formula, 



Example 20. Let A e 



IIAlb = 
l 4x3 be the matrix 



lim \\(AA>} 



k—>+oo 



( 3 
2 

-3 




2 
1 

25 
-1 



32^ 

36 
2 

1 ) 



Then, with the algorithm in IRemark 141 we get that the 2-norm is 48.46054603. Using the 
standard algorithm (the spectral radius of AA') we get the same number 48.46054603. 
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Example 21. An interesting example is the maximum of a bilinear form over S 1 x §°. Note 
that the domain is a cylinder of R 3 , so we can draw the whole situation. Take, for example, the 
bilinear form 

M 2 xl — > K, i{x,y) = 4x iy + 2x 2 y. 
The maximum of ( over S 1 x §° is the 2-norm of the vector (4, 2), that is, 

||(4,2)|| = V20 = 4.472135954. 

Let's compute this using the algorithm in IRemark 14\ First of all note that the gradient of ( 
determines a vector field over the cylinder, and the iteration follows the arrow. Over the ending 
point of the iteration, the flow is orthogonal to the surface. This means, that we have reached an 
extreme, 

max e(x,y) = 4.472135953. 

IWl=bl=i 



Let's give now some applications of the general algorithm, see Appendix A 



Remark 22. The first interesting application of the algorithm in Appendix A is to the theory of 
entanglement. It is of interest to find the maximum of the form (p, -} over the space of separable 
states. The matrix p is called a state if it is hermitian, p > and tr(p) = 1. It is easy to see that 
the space of states is a convex set and is generated by the matrices of the form pip\ where p,- is 
a column vector of norm one in a finite dimensional vector space 'H, p, e 'H, ||p,|| = 1. The 
general theory says that when we are working with two particles, we need to consider the space 
of states over the tensor product 'H = 'Hi ® 'Hi. In this situation a state is called separable if 
it is a convex combination of the form 2 fl;V; ® w,-, where v,- is a state of 'Hi and w; a state of 
7^2. Let's call SepCH) the convex space of separable states. It is true that the space of separable 
states is a convex set generated by the matrices of the form xx* ® yy* , where x € *H\, y e "Hi and 
IWI = \\y\\ = 1. Then 

max (p, -) = max (p, xx* ®yy'}. 

SepCH) IWHIvIN 

Note that the form is not bilinear in x nor in y. Rewriting the state p in the form p = 2 A-iPiPl 
with (pi,Pj) = 0, ||p,|| = 1, and using the equality (p;, x ®y) 2 = (piP, ,xx^ ® yy 1 ), we get 

max (p, -) = max || V V^KP;, * ® ^A'll 2 - 
SepCH) IWHIyIN ^ 

The resulting map 2 ypti(pi, x ® y)pt is bilinear in x and in y and our algorithm is capable to 
maximize it. See the next example. 

Example 23. Suppose that we are working with the following state p in R 4 = R 2 ® R 2 , 

0.242894940524649938 -0.123994312358229969 -0.0712215842649899789 0.219784373378769966 

-0.123994312358229969 0.0888784895376599772 0.111143109132249979 -0.0627261109839499926 

-0.0712215842649899789 0.111143109132249979 0.361255602168969903 0.0603142605185699871 

0.219784373378769966 -0.0627261109839499926 0.0603142605185699871 0.306970967813849916 

We choose to work over the real numbers to make the exposition clear, but all the results can be 
adapted to work with hermitian matrices instead of symmetric matrices. Using Cholesky and the 

Singular Value Decomposition Algorithm we have p = 2 ^iPipj, 

A { = 0.5435016101, A 2 = 0.4146107959, 
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A 3 = 0.04113792919, A A = 0.0007496649711, X\ + A 2 + A 3 + A 4 = 1. 



Pi 



P3 



-0.656481390369177854^ 
0.326643787198963642 
0.245965753592146592 

-0.633906040705653040 ) 

-0.444223726945872255^ 
0.546710519336902400 

-0.399219012690601172 
0.586853532298337144 



92 ■ 



P4 



-0.0253829550629408562^ 
-0.209402292907094082 
-0.881013881627254691 
-0.423463015737623016 ) 

-0.609 141 338368644^ 
-0.741998735885678107 
0.0627659808038323054 
0.272846362424434330 



With the algorithm in | Appendix A| to the trilinear form 

£(x,y,z) = J] 4Xi{p h x®y){ Pi ,z), (x,y,z) eS'xS'xS 3 



we get that the maximum is 0.7228016991. Finally 



max (p, -) = 

Sep(H) 



1 2 v^fo- 



)y) Pi \\ 2 = max \({x,y,z)\ 2 
IWI=llyll=llzll=i 



0.7228016991 z = 0.5224422962. 

Note that if p is separable then (p,p) < maxs ep (f{){p,-). This is not the case, but for example, 
the following state is not separable (it is called entangled state), 



0.168106937369559950 
-0.190509527669719958 
-0.200004375511779936 
-0.0690454833860399825 



-0.190509527669719958 
0.257651665981429912 
0.267759084652009926 

0.0985801483325399742 



-0.200004375511779936 
0.267759084652009926 
0.320790216378169901 
0.194053687463299957 



-0.0690454833860399825 
0.0985801483325399742 
0.194053687463299957 
0.253451180300149959 



We have <p,p> s 0.6620536187 ^ 0.4862909489 = msx SepW {p, ->. 

Remark 24. Our final application is the ability to find numerically the closest rank-one tensor 
of a given tensor. In the article [8], the authors considered the problem of finding the best rank- 
r approximation of a given tensor. They proved that for r > 1 the problem is ill-posed, but 
when r — 1 the problem has a solution, 18l 4.5], Here we find the solution. Let's prove that a 
computation of the absolute maximum of C over a product of spheres gives the closest rank-one 
multilinear form to t . A rank-one multilinear form is a product of linear forms, l\ . . . ( s , where 
(i : — » K, 1 < i < s. We choose to do this remark about multilinear forms, but dually, the 
same is true for tensors. 



For simplicity, we do the proof for a trilinear form. The proof is similar in the multilinear 
case. Consider the affine Segre map (it is not an isometry) 



r+ i x Rm+ i x w+ \ 



TJ/7+1 ,o, TTT>m+l 



(x,y,z) 

Using the usual inner product in the tensor product, we identify 



D« + l ,<-, TTT>m+l 



v+ l s (R «+1 ^ R m + l 



>S+l-jV 



x<8>y<8>z. 



x ®y ® z — > (x ®y ® z, -), (x ® y ® z, a ® b ® c) = (x, a)(y, b)(z, c). 
We can identify the following three different notations 



€{x,y,z) - t{x®y®z) - {(,x®y®z). 
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The first equality identifies a trilinear form with a linear map I : R" +1 ® M"' +1 ® R I+1 — > R. The 
second equality identifies, under the isometry (M" +1 ® M" ,+I ® R' s+1 ) v s R" +1 ® M" ,+1 ® R I+1 , the 
linear form € with the tensor t e R" +I ® R m+1 ® R s+1 . 

Let § be the immersion of §" x x §' 5 under the Segre map, 

S = {{x%y®z,-)\ \\x\\ = \\y\\ = Hzll = 1} c (M' ,+ 1 ® M"' +1 ® M' 5+I ) v . 

Then, for all <p — (x ® y ® z, -) € §, we have 

IK - 0II 2 = (IKII 2 + ll^ll 2 - 2((, (/>)) = IKH 2 + 1 - 2£(x,y,z). 

In other words, a local maximum of £ is a local minimum of the distance function, ||£ — 0||. 
Let B be the image, under the Segre map, of a product of balls, 

1 = {(x ® y ® z, -> : ||jc||, ||y||, ||z|| < 1} c (M' ,+ 1 ® R m+1 ® R I+1 ) V . 

Note that the elements of B are rank-one multilinear forms. It is easy to see that B is compact 
and convex, so the distance from I to B is achieved in § (the border). In other words, the closest 
rank-one multilinear form to I is an element of S. Summing up, a computation with the algorithm 
in lRemark 18l of the absolute maximum of t, gives a closest rank one multilinear form to t. 

Example 25. Let t : R 2 x M 2 x R 2 x M 2 — > R be a multilinear form, 

£(x,y,z,t) =Ax\y\zih +6x 2 yiZifi + xiy 2 ziti + lx2yiz\t\ -5x\y\Zit\- 
6x 2 y\Z2t\ - 5x\y 2 Z2t\ + 2x 1 y 1 zi? 2 - ~ix 2 y\Z\t 2 - lx x y 2 z\t 2 + 9x 2 y 2 z\t 2 - 
9x\y\z 2 t 2 - 9x 2 y l z 2 t 2 - 6x\y 2 z 2 t 2 + 8x 2 y 2 z 2 t 2 . 

Using the algorithm in lRemark 181 we get that the closest rank one multilinear form is 

(x(x)£ 2 (yYi(z)Ut), 

h(x) = 0.4799354720X! - 0.87730379 18x 2 

£ 2 (y) = 0.2732019392y!-0.9619567040y2 

e 3 (z) = 0.7563638894z! + 0.65415 11043Z2 

U(f) = 0.32609483 15f i +0.9453370622f 2 . 

The value of the absolute maximum of I is 16.71262553. 

Example 26. Let v € M 2 ® R 3 be the following tensor 

v = 4xi ® yi — 9x 2 ®yi + 2xi ® yi + x 2 ® y 2 — 5x\ ® y^ — lx 2 ®y3. 

Using the algorithm in lRemark 18l we get that the closest rank one tensor is 

(0.01162554952xi + 0.9999324213x 2 ) (-0.7821828869yi + 0.08939 19925 ly 2 - 0.6 166027924^) 

In this case, we can check this result. The first singular vectors of the matrix 

'4 -9^ 

2 1 , 
1-5 -V 

are 

(0.01162554952,0.99993242102), (-0.7821828866,0.08939199251,-0.6166027924). 
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A. General Algorithm for a multilinear form. 

Let's give an algorithm to find the maximum value of generic multilinear map over a product 
of spheres, 

I : K" 1 X ... X R"' — > K" r+1 , max ||£(*i,...,* r )||. 
Recall from lSection 2l that we may assume that i is a multilinear form, 

?: K' !l x ... x M" r x -> R, max . . . , x r+1 )|. 

I|a'iII=...=I|a-,. + ,||=1 

The following is a pseudocode in the trilinear case. We choose to work with Grobner Basis 
because it is implemented in most computer algebra systems (Maple, Macaulay2, Singular). In 
J2], the authors proposed an algorithm without the need of Grobner Basis. See also [9, 2.3.1]. 

Input: A generic trilinear form f:l' +l x8" +1 xr +l -» K. 

Output: The absolute maximum (x, y, z) e S" X §"' X S s . 

1. Compute the system of trilinear equations of [Proposition 16] 

2. Add the polynomial equations ||jt|| 2 = ||y|| 2 = ||z|| 2 = 1. 

3. Compute a Gbbner basis for the resulting system, I. 

4. Find a basis B of C[x\, . . . ,x n ,yi, . . . ,y m ,zi, ■ ■ ■ , Zj]/-T. 

5. Compute the multiplication matrix of t, 

6. Return the magnitude of the maximum real eigenvalue. 

Let's give an implementation of the algorithm in Maple. The code computes the maximum 
of a trilinear form over S" _1 x S" I_I x S . The reader may change the values of n, m and s and 
the trilinear form, to get different examples. 

> restart ; with (Groebner) : with(linalg) : 

> n:=2:m:=2:s:=2: 

> L : =6*x [1] *y [1] *z [1] +3*x [2] *y [1] *z [1] -6*x [1] *y [2] *z [1] +16*x [2] *y [2] *z [1] - 
14*x [1] *y [1] *z [2] -15*x [2] *y [1] *z [2] -ll*x [1] *y [2] *z [2] +8*x [2] *y [2] *z [2] ; 

> #Step 1 and 2 

> J:={add(x[i]-2,i=l. .n) -1 , add(y [j] ~2 , j=l . .m) -1 ,add(z [k] '1 ,k=l . .s)-l, 

> seq(seq(x[i]*diff (L,x[j])-x[j]*diff (L,x[i]) , j=l. . i-1) ,i=l. .n) , 

> seq(seq(y[i]*diff (L,y[j])-y[j]*diff (L,y[i]),j=l. .i-l),i=l. .m) , 

> seq(seq(z[i]*diff (L,z[j])-z[j]*diff (L,z[i]) , j=l. .i-1) ,i=l. .s)}: 

> #Step 3 

> G:=Basis(J, 'tord') : 

> #Step 4 

> ns ,rv : =NormalSet CG, tord): 

> #Step 5 

> mulMat : =evalm(evalf (Mult ipli cat ionMatrixCL,ns , rv,G, tord) ) ) : 

> #Step 6 

> max (op (map Cabs , map (Re , {eigenvalues (mulMat )>) ) ) ) ; 

The following is a table that shows the time, in seconds, used to compute the maximum. In 
the first column appears different values of (n, m, s), in the second, the time used to compute the 
Steps 1 through 4 and in the third, the total time of the algorithm. We run a Maple 1 1 session on 
a 2. 1GHz CPU, with 2GB of memory, 



(«, m, s) 


Steps 1-4 


Total time 


(2, 2, 2) 


0.03 


0.33 


(2,2,3) 


0.05 


0.79 


(2,2,4) 


0.09 


0.99 


(2,2,5) 


0.14 


1.20 


(2,3,3) 


0.31 


7.13 


(2,3,4) 


0.89 


30.03 


(3,3,3) 


5.06 


397.28 
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Note that the computation of the multiplication matrix using Grobner basis, requires most of the 
time. In J2l, the authors proposed an algorithm to compute the multiplication matrix of £ directly; 
without Grobner basis. 
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